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ABSTRACT 

Magnetic fields play an important role in star formation by regulating the removal of 
angular momentum from collapsing molecular cloud cores. Hall diffusion is known to 
be important to the magnetic field behaviour at many of the intermediate densities 
and field strengths encountered during the gravitational collapse of molecular cloud 
cores into protostars, and yet its role in the star formation process is not well-studied. 
We present a semianalytic self-similar model of the collapse of rotating isothermal 
molecular cloud cores with both Hall and ambipolar diffusion, and similarity solutions 
that demonstrate the profound influence of the Hall effect on the dynamics of collapse. 

The solutions show that the size and sign of the Hall parameter can change the size 
of the protostellar disc by up to an order of magnitude and the protostellar accretion 
rate by fifty per cent when the ratio of the Hall to ambipolar diffusivities is varied 
between —0.5 ^ Vh/va ^ 0.2. These changes depend upon the orientation of the 
magnetic field with respect to the axis of rotation and create a preferred handedness 
to the solutions that could be observed in protostellar cores using next-generation 
instruments such as ALMA. 

Hall diffusion also determines the strength and position of the shocks that bound 
the pseudo and rotationally-supported discs, and can introduce subshocks that further 
slow accretion onto the protostar. In cores that are not initially rotating (not examined 
here), Hall diffusion can even induce rotation, which could give rise to disc formation 
and resolve the magnetic braking catastrophe. The Hall effect clearly influences the 
dynamics of gravitational collapse and its role in controlling the magnetic braking and 
radial diffusion of the field merits further exploration in numerical simulations of star 
formation. 
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1 INTRODUCTION 

Low-mass stars form by the gravitational collapse of molec- 
ular cloud cores over many orders of magnitude in size 
and density. Cores form within molecular clouds as a re- 
sult of turbulent fluctuations, and they gradually contract 
as ambipolar diffusion erodes the magn etic support (e.g. 
iTasker fc T"an.,200a : Federrath et al.ll201(]| '). Turbulence may 
also support the core against col lapse and its decay can aid 
in triggering star formation (e. g. iMac Low fc Klessenll20o3 : 
iBallesteros-Paredes et al]l2007l 'l. The core becomes unstable 
when the mass-to-flux ratio exceeds the critical value 



(1) 



where is a dimensionless (in cgs units) numerical coeffi- 
cient that depends upon the magnetic field and density dis- 
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tributions llMestel fc Spitzerl Il956l : iMouschovias fc Spitzerl 
ll976l : lMcKee et al.lll993l )" 

The core collapses dynamically into what is termed a 
"pseudodisc" (Galli fc Shu 1993a, b), which has a flattened 
shape due to the material falling in preferentially along the 
magnetic field lines. The field is effectively frozen into the 
cloud, and the infalling material is defiected by the field 
lines towards the equatorial plane. The pseudodisc contracts 
dynamically in the radial direction, dragging t he field lines 
into a split monopole (hourglass) configuration (iGalli fc Shul 
that is consistent with observat ions of the polarisa- 



mtoa spl 
Il993al lbl) 



tion of dust continuum emission in cores |Cortes fc Crutcher 
20061 : iGirart et aLll2006l : IConcalves et al.ll2008l : lAttard et al 



20091 ). The build up in magnetic pressure acts as an imped- 



iment to further collapse, however, the magnetic tension in 
the envelope never suffices to suspen d the envelo pe against 
the gravity of the growing protostar (|Allen et al. 2003a b). 

Wit hin the pseudodisc ambipolar diffusion becomes im- 
portant (jPesch fc Mouschoviad boOlh . leading to a decou- 
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pling of the magnetic field from the neutral gas, which 
takes place within an outwardly-propagating MHD shock 
(jLi fc McKeelll996h . This shock (referred to as the "mag- 
netic diffusion shock" in this work) i s a continuous transition 
in the magnetic field and density l|Krasnopolskv fc Konigll 
12002. hereafter KK02f ). inwards of which the neutral ma- 
terial falls in at a near-free fall speed due to the reduced 
magnetic support. The centrifugal force becomes important 
and triggers the formation of a hydrodynamic shock that 
strongly decelerates t he infalline mat ter and allows a Ke- 
plerian disc to form l|Shu et al.l ll987l ). A disc wind or jet 
may form, launched from the inner regions of t he collapse 
l|K6nigllll989l : lTomisakall2002l : lAUen et al.ll2003bl ). 

Angular momentum is removed from the pseudodisc by 
the twisting of magnetic field lines, which transport angular 
moment um from the inner parts of th e core towards its outer 
regions (jBasu fc Mouschoviad [l994l ') . The amount of mag- 
netic braking affecting the collapse, and hence the existence 
and sign of the Keplerian protostellar disc, is determined by 
the coupling of the field to the charged particles and the drift 
of these against the neutral component in response to the 
electric field in the neutral rest frame. The Lorentz force 
is transmitted to the neutral gas through the drag forces 
caused b y collisions between the n eutral and charged parti- 
cles (e.g. lKonigl fc Salmeronll201ll ). 

Simulations of star formation typically approximate the 
magnetic field behaviour by ideal magnetohydrodynamics 
(IMHD), where the mass-to-flux ratio is held constant and 
the magnetic field is regarded as being frozen into the 
neutral medium (e.g . iGalli et al.l l2006l : iMellon fc lH l2008l : 
iMachida et al.ir2008al V In this situation the magnetic field 
and the particles move together in the collapsing flow, how- 
ever this simpliflcation only truly applies in the outermost 
regions of gravitational collapse where the density is low. If 
IMHD were to hold true throughout the collapse the mag- 
netic flux in the star would be lO'^-lO^ times larger than that 
observed in young stars (this is the "magnet ic flux problem" , 
described in IChandrasekhar fc Fermill 19531) . 

As the density in the core increases flux freezing breaks 
down, and the relative drifts of different charged species with 
respect to the neutral particles delineate three magnetic dif- 
fusivity regimes: 

• the Ohmic (resistive) diffusion limit, which dominates 
in high density regions where the ionisation fraction is low. 
The ions and the electrons frequently collide with the neu- 
trals over the electron gyration period, and the magnetic 
field is decoupled from all charged particles. Ohmic diffu- 
sion is important in the innermost regions of the protostel- 
lar disc where the density and collision al rates are high (e.g. 
IShu et al.llioo e': Machi da et al.ll2008bl ). 

• the ambipolar diffusion limit, which dominates in re- 
gions of relatively low density where the fractional ionisa- 
tion is high, causing the ionised component to drift with the 
field through the neu trals. Ambipo lar diffusion is dominant 
in molecular clouds (|Wardlell2007l '), in protostellar discs at 
radial distances beyond ~ 10 au an d close to the s urface of 
these discs nearer to the protostar (|Salmeronll2009l ) . 

• the Hall diffusion limit, which dominates in the inter- 
mediate regimes between ambipolar and Ohmic diffusion. 
The more massive particles such as ions and charged dust 
grains are decoupled from the magnetic field and are instead 



coUisionally-coupled to the neutral gas. Hall diffusion is ex- 
pected to dominate in many regions of molecular c louds as 
they undergo gra vitational collapse (|Wardlell2"004al ). and in 
protostellar discs (|Sano fc Stone|[2002a 



In the Hall limit the magnetic response of the disc is 
not invariant under a global reversal of the magnetic field 
l| Wardle fc Nd 1 19991 ) as Hall diffusion twists the field and 
changes the angular momentum of the neutral fluid; in a 
rotationally-supported disc this causes the gas to fall in- 
wards if it loses angular momentum, or outwards if it gains 
it. Ambipolar and Ohmic diffusion, however, always cause 
the fleld to move in the radial direction against the flow of 
the neutrals - reversing the direction of the field does not 
affect the direction of the field diffusion. 

The actual magnitude and type of coupling that occurs 
between the fiuid and magnetic field in molecular clouds 
and protostellar discs is uncertain due to the difficulty in ob- 
taining detailed observations in these regions, particularly of 
the magnetic field . Calculations o f the ionisation equilibrium 
and resistivity bv lWardl3 (|2004al ) suggested that Hall diffu- 
sion is important and may dominate the magnetic field be- 
haviour at many of the densities and field strengths encoun- 
tered in molecular clouds and protostellar discs. In partic- 
ular, the Hall term is significant for molecular gas densities 
in the range ~ 10^-10^^ cm~^ (when B scales as _B oc n^^'^), 
although the presence and distribu tion of grains compl icates 
the calculation of the diffusivities (| Wardle fc NdlT999l ). 

In simulations of star formation where the magnetic 
field behaviour is not governed by I MHD it is usually am- 
bipolar diffusion tha t is included (e .g. Ciplek fc Konigl [l998l : 
lAdams fc Shul l2007l : iMellon fc lH |2009| ). Ohmic diffusion 



is important in the innermost regions of the collapse as 
the^ density builds up, particularly when p > 10^^ cm~'^ 
(iShu et al..2006 ). It seems likely that star formation requires 
all of these processes to some degree, however it is only re- 
cently that simulations have be en performed w ith more than 
one of these processes included (|Li et al.ll20Tll ). and Hall dif- 
fusion is almost always overlooked. The nature of the cou- 
pling determines the magnetic field direction as the field 
lines emerge from the surface of the protostellar disc, which 
in turn controls the amount of material that is able to slide 
along the field lines and be flung outwards from the surface 
in disc -driven wind models (| Wardle fc Koniglll 19931 : IWardld 
l2004bl ). 

In this paper we construct a semianalytic model of grav- 
itational collapse in order to demonstrate the importance 
of Hall diffusion and its influence on the magnetic braking 
catastrophe, in which the magnetic braking affecting a col- 
lapsing flow is such that all angular momentum is removed 
from the flow and no Keplerian disc may form. The paper 
is organised as follows: in §2 we describe the formulation of 
the self-similar collapse equations, the assumptions govern- 
ing this approximation and the boundary conditions of the 
core; the numerical procedures are outlined in §3; in §4 we 
present the similarity solutions; the effect of the Hall term 
on the solutions and the magnetic braking catastrophe is 
discussed in 55 and our conclusions are stated in 56. 
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2 FORMULATION 

The goal of this work is to construct a semianalytic model 
of gravitational collapse similar to that of KK02, including 
terms for Hall diffusion in the equations for the magnetic 
field diffusion and braking. This allows the calculation of 
similarity solutions that show the importance of Ifall diffu- 
sion in molecular cloud cores and collapsing flows, as well as 
comparisons between the influence of the Hall and ambipolar 
diffusion terms. Following KK02, the magnetohydrodynamic 
equations for the isothermal system are given by 



|+V.(pV) = 0, 



(2) 



P^+p(V- V)V = -VP + pg + J X B, (3) 



V-B = 0, 



(4) 
(5) 



and 



f = Vx(VxB) 

-Vx [r7(VxB) +?7h(VxB)xB + r7A(VxB)J, (6) 

where p is the gas density, V the velocity fleld, P the pres- 
sure, g the gravitational fleld, $ the gravitational potential, 
c the speed of light, J the current density, B the magnetic 
field, and rj and TjH,A are the diffusion coefficients for the 
Ohmic, Hall and ambipolar terms in the induction equation. 

We assume that the collapse is axisymmetric to simplify 
the calculations, and the magnetic field is taken to be aligned 
with the axis of rotation. This is in contrast with some obser- 
vations of protostellar systems such as the binary NGC 1333 
IRAS 4A where the axis normal to the binary envelope lies 
betw e en the outflow and magnetic fleld axes (iGirart et al.l 
I2OO6I : lAttard" et ahlfioogl ). however most observations seem 
to show that the magnetic fleld takes on an hour glass shape 
that i s aligned with the axis of rotation (e.g. IVink et al.l 
I2OO5I : ICortes fc Crutcheil [ 20061 ) . Alignment is required by 
the assumption of axisymmetry, as misalignment introduces 
three-dimensional effects that cannot be modelled here. 

Using cylindrical coordinates, the mass, radial momen- 
tum and angular momentum conservation equations, as well 
as the hydrostatic equilibrium equation, the solenoidal con- 
dition and the vertical component of the induction equation 
are, under the assumptions of isothermality (that P 
and axisymmetry: 
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where Qr and g, are the radial and vertical components 
of the gravitational field, and Cs is the isothermal sound 
speed given by Cs = {ksT /rUnY^^ ~ 0.19 km s~^ (with ks 
the Boltzmann constant, T the gas temperature, typically 
taken to be 10 K, and m„ the mean mass of a gas particle). 
The assumption of isothermality breaks down due to radia- 
tive trapping w hen the central density reaches ~ 10^" cm""^ 
l|Gaustadlll963 ). which occurs on scales r < 5 au for a typi- 
cal simulation. It is expected that isothermality shall break 
down in the innermost regions of our solutions, however, as 
thermal stresses do not play a significant role in the larger- 
scale dynamics isothermality is not expected to introduce 
large errors into the calculations. 

The disc is assumed to be thin b ased upon the results of 
Mouschovias and collaborators (e.g. iFiedler fc Mouschoviad 
1I992 , ,1993), which have shown that an initially-uniform, 
self-gravitating, magnetised molecular cloud core rapidly 
collapses along the magnetic field lines. This assumption al- 
lows us to further reduce the dimensionality of the problem 
by vertically-averaging the variables over the scale height of 
the disc, although it implies that processes that depend on 
variations in the density or the magnetic field with height 
within the disc cannot be included in the collapse calcula- 
tions. Effects such as turbulence or the interaction between 
active and dead zones in the disc are not expected to have 
a large effect on the overall dynamics of early collapse (al- 
though they are known to become important in some regions 
of protostellar discs once the adiabatic core and protostar 
have formed) and their exclusion is necessary to the self- 
similarity of the solutions. 

The vertical averaging is performed as in KK02. As only 
the vertical component of the induction equation and our 
equation for the magnetic braking differ from theirs, we per- 
form the integration of these in Appendix |X] and refer the 
reader to their appendix A for the others. The quantities 77, 
r)H / B and tia/B'^ are approximated as being constant with 
height, as are the radial velocity, the azimuthal velocity and 
the radial component of gravity. 

The disc is threaded by an open magnetic field pos- 
sessing an even symmetry, so that Br = B^ = at the 
disc midplane. The radial and toroidal field components are 
taken to scale with height as 



Br{r,z) =Br,s{r) 
B^{r,z) = B^,s{r) 



H{r) 
z 

mr) 



(13) 
(14) 



where H{r) is the scale height of the disc and Br.s and B, 
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are the magnetic field components at the surface of the disc; 
these scahngs are motivated by the field configuration of a 
rotationally-support ed thin disc in which t he field is well- 
coupled to the gas jWardle fc Konig]|l993l '). As ambipolar 
diffusion and Hall diffusion become more important in the 
inner regions of the disc where the field is less well-coupled to 
the gcis, this approximation is no longer adequate, however 
none of the dominant terms in the equation set depend upon 
the particulars of the vertical variation of the field within the 
disc (KK02), so it remains reasonable to adopt these scalings 
across the domain of self-similar collapse. 

The surface density of the pseudodisc is defined by the 
expression 



r 



pdz = 



(15) 



assuming that the density is constant with height within 
the disc, and the specific angular momentum is defined by 
J = rV^. We neglect any mass loss due to a disc wind. 
The radial components of gravity and the magnetic field are 
calculated using the monopole expressions 



and 



9r 



GM(r) 



27rr2 ' 



(16) 



(17) 



where the enclosed mass M{r) « when the central mass 
dominates and ^ is the magnetic fiux enclosed within the 
radius r. These simplific ations were also used by KK02, as 
IContopoulos et all (| 19981 ') found that these expressions give 
values of the gravitational force and magnetic field that are 
near enough to those found using an iterative method that 
they do not introduce significant errors into the calculation. 

The vertical angular momentum transport above and 
within the disc is achieved by magnetic braking, especially 
during the dynamic collapse phase inwards of the magnetic 
diffusion shock and in the rotationally-supported inner disc. 
The approach to modelling the magnetic braking used in 
this work, which is desc ribed in more detail in Appe ndix Rl 
is adapted from that of lBasu fc Mouschoviaj l| 19941 ) for the 
pre-point mass formation collapse phase. This formulation 
is not well-defined in the innermost rotationally-supported 
regions of the pseudodisc, where the calculated magnetic 
braking becomes stronger than is expected and the angular 
momentum transport is expected to be dominated by a disc 
wind (which is not included but discussed further in ij5]). 

A cap is then placed upon the azimuthal magnetic field 
component in order to ensure that it does not greatly exceed 
the vertical component; balancing the torques on the disc 
then allows us to define 



Brh, = — min 



TTr'^ V VA.cxt 



(18) 



where 5 is a constant parameter {— 1 in our solutions) that 
limits the magnetic braking, Q is the local angular veloc- 
ity, and fib is the background angular velocity of the cloud, 
which is small compared to that of the core. VA.cxt is the 
Alfven wave speed in the external medium, which is param- 
eterised with respect to the sound speed using 



where again a is a constant parameter of the model. We 
take a = 0.08 in the solutions presented here, which is a 
reason able approximation as the observations of ICrutcheJ 
(Il999h indicated that VA,cxt ~ 1 km s over many orders of 
magnitude in several molecular clouds. The negative sign of 
_B^_s and positive sign of S are based on the assumption that 
there is no counter-rotation of material in the collapse; were 
the fiuid counter-rotating then the signs of 5 and B^^s would 
both change. The azimuthal drift velocity of the magnetic 
field is averaged over the disc scale height in Appendix 1X1 to 
give the final form of Equation IA37I shown below. 

The disc equations are further simplified by recognising 
that the thin disc approximation implies that terms of order 
0{H/r) are small in comparison to the other terms and 
can then be dropped from the equations. As in KK02, the 
only term of order 0(H/r) that is kept is the combination 
[_B,.,s — H{dBz/dz)], which occurs in the radial momentum 
equation and is important in refining the structure of the 
magnetic diffusion shock. This term is then retained in all 
of the equations in which it appears. 

Taking all of these into account then gives the simplified 
set of equations: 



dVr , dVr 



+ 



S dr 2nT. 



+ 



dt dr 27rE 



(20) 



(21) 



(22) 



= -GS' + ^^ + ^(B.% + BL), (23) 



2H 2 



E.d^ = .rHKB. - - If B.B,,. - '-^B^.^bI 

(24) 



and 



B^^s = - min 



'^a r\p Bz 



\5Bz 



(25) 



s/Va.c 



(19) 



along with Equations 1161 and 1171 
2.1 Self-Similarity 

At any instant in time the collapse solutions look like 
stretched versions of themselves at previous times; this 
fractal-like behaviour is referred to as self-similarity. The 
pseudodisc forms as a collapse wave (referred to as the mag- 
netic diffusion shock) propagates outwards at the speed of 
sound. The self-similarity of the waves of infall occurs be- 
cause of the lack of characteristic time and length scales in 
the flow. 

Gravitational collapse occurs over many orders of mag- 
nitude in radius and density, so that the point mass has 
negligible dimensions in comparison with the accretion flow. 
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The only dimensional quantities that effect the flow are the 
magnetic field B, the diffusion coefficients 77 and »?h,a, the 
gravitational constant G, the isothermal sound speed Cs , the 
local radius r and the instantaneous time t; this means that, 
except for scaling factors, all of the flow variables may be 
written as functions of a similarity variable deflned by 



r 

Cst' 



(26) 



KK02 noted that for a typical value of the sound speed (cs = 
0.19 km s~^ at T = 10 K), x — 1 corresponds to a distance 
of r 6 X 10^^ cm (400 au) when t = lO" yr (which is the 
characteristic age of a Class YSO) and to a distance of 
r ^ 6 X 10^^ cm (4,000 au) when t = 10^ yr (the age of a 
Class 1 YSO). The Class YSO IRAM 04191 has a dense 
inner disc-like structure that re sembles a tilted ring with an 
average radius of ro ~ 1400 au (|Lee et al.lboosi ') — this is of 
the same order of magnitude as the centrifugal shock radius 
in the disc-forming solutions at the same age. 

The physical quantities are expressed as the product of 
a nondimensional flow variable that depends only upon x 
and a dimensional part constructed from Cs, G and t: 



S(r,t) 

gr{r, t) 

Vr{r,t) 
H{r,t) 

J{r,t) 
M{r,t) 

M{r,t) 
B(r,t) 
*(r,t) 
and riH,A 



\2-KGt 



CsU{x), 

Csth{x), 

CsV{x), 



m{x), 



rn{x), 



(g^/h) 

V GV2" 

2, / 
CstVH,A- 



h{x), 



(27) 

(28) 

(29) 
(30) 
(31) 
(32) 

(33) 

(34) 
(35) 
(36) 
(37) 



These equations have the same form and use the same nota- 
tion as those in KK02, with the addition of extra diffusion 
coefficients to model the magnetic field more completely. 

The Ohmic and ambipolar diffusion terms scale to- 
gether, to a zeroth-order approximation, as they possess a 
similar dependence upon B and appear in the induction 
equation and the equation for the azimuthal field compo- 
nent multiplied by the same magnetic field terms. Because 
the field within the thin disc is effectively vertical, both am- 
bipolar and Ohmic diffusion influence the fleld drift in the 
same manner. While one type of diffusion may dominate 
over the other at any individual point in the disc (in general, 
ambipolar diffusion in the outer regions where the density is 
low and Oh mic diffusion in the inner regions where the den- 
sity is hig h: IWardlell2007l l. only one term is needed in order 
to study the change in the disc behaviour introduced by the 
Hall diffusion term that is of most interest in this work. The 
Ohmic and ambipolar diffusion terms are combined into a 
single term parameterised by the dimensionless constant fjA , 
referred to as the ambipolar diffusion parameter. 



The ambipolar diffusion coefficient in a molecular cloud 
core without grains is given by the equation 



47r7pip 



(38) 



where 1/ypi = T„i is the neutral-ion momentum exchange 
timescale, parameterised as 



VA 



(39) 



the nondimensional ambipolar diffusion parameter fjA is a 
constant of the model (simply denoted rj in KK02). tja/B^ 
is then self-similarised using the scalings above to give 



62 "'''*cr3/2' 



Va 



(40) 



it is important to note that the self-similarity of the solution 
depends upon the relationship pi oc p^''^ ■ 

For grains with radius a = 0.1 pm in a cloud where the 
temperature is 10 K and the cosmic ray ionisation rate is 
^ = 10~^^ s~^, simulations typically assume the ion den- 
scales as Pi oc pl/^ when 10* < nn 5- 10^ cm~^ 



sity ^ _ 

(|Elmegreenlll979l: iKamava fc Nishill2000l). This is an over- 
simplification, as ICiolek fc MouschoviasI |l998l ) showed that 
for typical cloud and grain parameters the proportionality 
of the ion density cannot be parameterised by a single power 
law exponent, but it is still a reasonable and widely-adopted 
approximation to the ion density in collapsing cores on 
scales > 10^ au (see e g.lShu et al.(ll987l: ICalh fc Shulll993al : 
ICiolek fc K5nig||l998l : ]ContoDoulos et al.lll998l . KK021 . 

As a matter of pragmatism, a similar scaling with re- 
spect to the density and scale height is adopted for the Hall 
diffusion parameter, rjH- By stating that the self-similar Hall 
diffusion coefficient scales as 



(41) 



where f/n is the constant nondimensional Hall diffusion pa- 
rameter used to characterise the solutions, the ratio of the 
nondimensional ambipolar and Hall diffusion parameters be- 
comes the most important factor in determining the mag- 
netic behaviour of the similarity solutions. In truth, the Hall 
diffusion coefficient could be scaled with respect to the den- 
sity and field strength by multiplying the nondimensional 
Hall parameter by any function of the similarity variable x 
and the fiuid variables. This topic is discussed in more detail 
in ij5l where an alternate scaling is proposed for future work 
on the self-similar collapse model. The scaling of r/'u given in 
Equation [41] is appropriate for a molecular cloud core with 
grains acting as the dominant positive charge characters. 

For convenience the variable w = x — u is used to sim- 
plify the equation set. The similarity variables are then used 
to rewrite Equations 120)4241 in self-similar form: 



dx 
dm 
dx 
Ida 
a dx 

dj 



dx 



= xb. 



= 9+- 
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dz 
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+ 



xb,b, 



2 "(^,3 



(42) 
(43) 
(44) 
(45) 
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^-br,s^]h' + +bl^+a')h-2a = 0, (46) 



and 



ip — xwbz + fjHxb^^sbzbh 

, ~ i2,l/2 -3/2 

+ rjAxb^h ' a ' 



1/2 -3/2 



I— 

dx 



= 0. 



(47) 



These equations are augmented by the self-similar defini- 
tions 



(48) 
(49) 

(50) 

while the other magnetic field components are given by 



m — xwa, 
m — —xua, 

and g = -; 

X'^ 



br,i 



(51) 



and 



= — min 



X V dx 



1 + 



2afiAh}^'^ipb:, 



„2n-3/2 



■,5b. 



(52) 



These equations completely describe the collapse of the 
molecular cloud core into a pseudodisc and the accre- 
tion onto the central point mass (potentially through a 
rotationally-supported disc). These equations are the same 
as equations 20-32 of KK02 in the limit of fjn ~ 0, which 
allows direct comparisons to be made between the similarity 
solutions of both models. 



2.2 Outer boundary conditions 

The outer regions of the collapse are modelled by a set of 
power law relations in the similarity variable that describe a 
molecular cloud core contracting quasistatically under am- 
bipolar diffusion until it has just beco me supercritical and 
a poi nt mass forms at the centre (e.g. IShulll977l : IShu et al] 
Il987h . 

The definition of the similarity variable (Equation I26p 
means that the limit a; — >■ oo corresponds both to the outer 
edge of the core at r — >■ oo and the initial conditions of the 
collapse as t — ^ 0. We then describe the core as a singular 
isothermal sphere, which has the density profile p ex (e.g . 
iLarsonl 119691 : |Penston|[l969l : IWhitworth fc Summers! [l985h . 
so that the surface density is: 



T,{t = 0) = 



(53) 



27rGr 

where A is a constant determined by the initial accretion 
rate of the core. The infall velocity and accretion rate onto 
the core are constant and given by 



Vr{t = 0) = UoCs 

and M{t = 0) = ~^^. 



(54) 
(55) 



The numerical results o f lCiolek fc Konigll l| 19981 ) for the col- 
lapse of a rotating magnetic core with ambipolar diffusion 
showed that the accretion rate at point mass formation was 



M ~ 5 Mq Myr ^, corresponding to a nondimensional pa- 
rameter A ~ 3, which also matches to obse rvations of man y 
cores that show M G [1,10] Mq Myr"^ l|Lee et all boOll ). 
As the isothermal sound speed of a core at T = 10 K is 
equal to Cs = 0.19 km s~^, the nondimensional initial infall 
speed must be of order unit y to match obs ervations show- 
ing Vr = 0.05-0.10 km s"^ (|Lee et al.|[200ll ): to match the 
parameter of KK02 we adopt uq = —1. 

The rotational velocity of the initial core is spatially- 
uniform and given by 

Vcf, = vaCs (56) 

where vq is the dimensionless rotational velocity, which can 
be approximated by 

^ A^lbCs (57) 

using the depende nce of the core surface density 

and magnetic field as in iBasul 11993). The uniform back- 
grou nd angular velocity i s typically fli, ~ 2 x 1 0~^^ rad 
s-i (|Goodman et al.lll993l : iKane fc Clemens! 1 19971 ) and the 
background ma gnetic field is taken to be Brcf = 30 
(|Crutcher!ll999D . which gives vo = 0.15. This v alue is a factor 
of ten larger than that obtained in lBasu!(|l997l ), ho wever, the 
range of observed core velocities is vo G [0.01, 1.0] (|Lee et al.l 
2001). KK02 showed that the size of the inner Keplerian disc 
directly corresponds to the initial rotational velocity at the 
core; we take vq = 0.73 to facilitate comparison with their 
fiducial ambipolar diffusion solution. 

Finally, in the outer regions IMHD holds true, and the 
mass-to-flux ratio in the gas is constant. The core is just 
supercritical at 

where fiQ is the dimensionless mass-to-flux ratio, parame- 
terised with respect to the critical v alue for support against 
gravity (jNakano fc Nakamura![l97^ ). A value of ~ 2.9 is 
adopted here for compatibility with KK0 2; this value was 
obtained from the numerical simulations o f lCiolek fc Konigl! 
( 1998. ) and matches observations showing that cores are 
typicall y more than twice supercritical (Crutcher 1999; bu t 
see also ICrutcher et~all |2009( : iMouschovias fc Tassid 12009! ) . 
Equation |42] can then by used to show that 

2-^- (59) 



Mo 



this is equivalent to Equation 1581 so that only one of these 
may be used to calculate the similarity solution. 

In self-similar form the conditions at the outer boundary 
Xout take the form 



m 


— Axout , 


(60) 


a 


A 

:^out 


(61) 


b. 


a 

MO 


(62) 


and V 


= Uq. 


(63) 



2.3 Inner boundary conditions 

As the variables are integrated inward from the centrifu- 
gal shock they tend towards an inner asymptotic set of 
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equations describing a Keplerian disc around the protostar. 
These relations are found analytically by assuming that the 
variables take the form of power laws in x and solving Equa- 
tions |42}[52] for the exponents as x — >■ 0, under the assump- 
tion that there is no counter-rotation. In nondimensional 
form this disc is described by: 



m 



■ ai X 



-3/2 



2 



-3/2 




1/2 



^3/2 



br- 
and bd, 



= -5 b, 



(64) 
(65) 

(66) 

(67) 
(68) 

(69) 
(70) 
(71) 

(72) 

(73) 
(74) 



The diffusion constant / is a function of the ambipolar and 
Hall diffusion parameters, and is given by the equation 



(75) 



this definition shows how the Hall term is able to counter- 
act the ambipolar diffusion term in determining the surface 
density of the disc and the accretion rate onto the central 
protostar when the nondimensional Hall parameter fjH is 
positive, and add to the ambipolar diffusion if the Hall pa- 
rameter is negative. The characteristic diffusion parameter 
of the disc, /, must be positive in order to form a Keple- 
rian disc in solutions without counter-rotatiorQ; this places 
limits on the relative sizes of the diffusion parameters. 

These relations are derived elsewhere (Braiding & War- 
die, in prep): to summarise, the inner accretion disc is in Ke- 
plerian rotation in the same direction as the initial core, with 
the centrifugal force balancing the inwards pull of gravity. 
The accretion rate onto the protostar is constant and low, 
and accretion through the disc is determined by the total 
amount of magnetic diffusion, which removes radial support 
by the magnetic field and allows the gas to fall inwards. 

A second asymptotic solution also exists (KK02; Braid- 
ing & Wardle, in prep), in which the magnetic braking is so 
strong that most of the angular momentum is removed from 
the gas, which then free falls onto the protostar without 
forming a disc (the magnetic braking catastrophe, described 
in Similarity solutions matching onto this solution lie 
beyond the scope of the present paper. 



^ As discussed in §5, the Hall effect can cause counter-rotation of 
the disc; however, in all of our solutions the inner Keplerian disc 
is rotating in the same direction as the initial core. 



3 NUMERICAL METHODS 

The calculation of the similarity solutions is quite com- 
plex, as many of the derivatives are large, and sonic points 
and shock fronts must be calculated explicitly as they are 
encountered by the integration routine. Furthermore, sub- 
shocks can occur downstream of the magnetic diffusion and 
centrifugal shocks. Such a subshock occurred downstr eam o f 
the magnetic diffusion shock in the solutions of Q (|l998h . 
and although no subshocks appeared in the published solu- 
tions of KK02 they were observed in the unpublished col- 
lapse solutions discussed in that paper. Locating and inte- 
grating through these shocks requires careful monitoring of 
the integration and automatic intervention where necessary. 
The full details of the methods used t o compute the s imilar- 
ity solutions are described in §5.1 of iBraidind (|201ll ): only 
an outline is provided in this work. 

The problem is recast as a two-point boundary value 
problem in which the variables are integrated from a match- 
ing point Xm to both the inner and outer boundaries. We em- 
ploy a "shooting" routine, which modifies the values of the 
variables at Xm in order to zero the discrepancies between 
the integrated variables and their expected asymptotic val- 
ues at the boundaries. The matching point is located at a 
position Xc < Xm < Xd, where Xd and Xc are the positions of 
the magnetic diffusion and centrifugal shocks respectively, 
and is typically chosen just downstream of the magnetic dif- 
fusion shock at x,n ~ 0.3. 

The values of the variables m, a, j, tp and b, at Xm and 
the nondimensional accretion rate onto the central mass rric 
are initially unknown, and a poor guess of these can cause 
the integration (and by extension convergence on the true 
solution) to fail. An acceptable initial guess of the variables 
is obtained by performing a simplified integration from the 
outer boundary to the matching point. 



3.1 Initial guess at Xm 

The initial guess of the variables at the matching point is 
estimated by performing a simplified integration in which 
the induction equation has been replaced by an algebraic 
expression for the vertical magnetic field component. We 
follow the derivation of KK02 here, and in the ambipolar 
diffusion {fjH = 0) limit our equations reduce to theirs. 

The initial guess of nic is estimated from the mass 
plateau in the region of the magnetic diffusion shock (see 
Q. The outer edge of this plateau occurs at Xpi « |wo| 
(which is something of an overestimate as \u\ is typically 
larger than x here). The outer asymptotic behaviour is only 
just starting to break down in this region, and so the asymp- 
totic relations from H2.2l are substituted into Equation 1481 



■nipi ^ 2|«o|A. 



(76) 



The calculation of the similarity solutions is less sensitive to 
our estimate of rric than to the values of the other variables 
at the matching point. Equation [76] is an acceptable first 
guess of the accretion rate onto the central star. 

Save for during the various shock transitions where b, 
changes rapidly, the inequality br.s ^ h{db,/dx) holds true 
everywhere in the collapsing flow; as this term is always 
small compared to the other terms in the induction equation 
it may be dropped. The magnitude of the magnetic field b 
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is always of order ~ &r,s, so that it may be estimated by 
b ~ \/2b2 . The induction equation may then be written as a 
quadratic in 6^: 

^^1/2^-3/2 (^fj^^b^^^ + riAbr,s) bl - xwb, = 0. (77) 

The well-separated roots provide approximations to the be- 
haviour of bz on either side of the magnetic diffusion shock, 
with the prescription 



a 
MO 



(78) 



applying in the large x regime where flux freezing still mostly 
holds true and the mass-to-flux ratio is given by its initial 
value fj, — fiQ. This is equivalent to the initial condition 
for the vertical field component and although IMHD breaks 
down before the magnetic diffusion shock Equation [TS] re- 
mains a good approximation to the field in this region. 

The larger root gives the value of the vertical field com- 
ponent in the magnetic diffusion regime where x is small, ft 
is approximated by dropping the constant term in Equation 
[77] and solving for 6z to obtain 



high ~ — (^-^^ (^V^f]Hb4,,s + f]A 



(79) 



which reduces to equation 50 of KK02 in the fjH ~ limit. 
The thickness of the disc in this region is controlled by the 
magnetic squeezing, so that 

h^^. (80) 

The transition between the two approximations to fez 
occurs at the magnetic diffusion shock, Xd, which represents 
a continuous increase in the magnetic field strength, and 
though the matter is slowed in the post-shock region, w is 
near-constant and w > 1 throughout the magnetic diffusion 
shock itself. The magnetic field lines are compressed by the 
shock, resulting in an increase in the vertical and azimuthal 
field as the field lines are twisted up by the slowing of the 
compressed gas. Conservation of flux ensures that the radial 
field component does not change in the shock. On the up- 
stream side of the magnetic diffusion shock b^^s ~ —Sb^^iow 
(as can be seen in the solutions in Q , however downstream 
of the shock the azimuthal magnetic field component is given 
by fe</),s ~ —wSb^^iow In the region of the magnetic diffusion 
shock fe<^,s is then simplified into 



-Xd5br,s 



(81) 



as Xd < 1, fe0,s < fer,s when S — 1, justifying our omission of 
the azimuthal magnetic field component from the approxi- 
mations to the scale height and magnetic field amplitude. 

Then using these relations and the approximation that 
IMHD holds true in this region. Equation 1791 becomes 



fez, high ~ — (V2fiA — 25Xdf}H 
X ^ 



(82) 



applying the approximation that b^^high = fer.s gives an es- 
timation of the magnetic diffusion shock position: 



f]A ( IJ.0 , 



(83) 



which is equivalent to KK02's equation 58 in the rjH = Q 
limit. This value of Xd is typically accurate to 20 per cent of 



the true value, which is acceptable for estimating the vari- 
ables at Xm- 

In order to obtain the initial guess of the variables at 
Xm we integrate Equations I42H46I inwards from the outer 
boundary, using Equation [78] to approximate fez when x > 
Xd, and Equation 1791 for Xm < x < Xd- These are used to 
integrate the full set of equations in both directions from Xm', 
the deviation of the results from the expected asymptotic 
boundary conditions is minimized by the shooting routine, 
which tweaks the variables at Xm to reduce the discrepancies. 
Once the initial values at close to the true solution 

the shooting routine will converge quadratically. 

Integration of the full equation set in the outwards di- 
rection is usually performed without difficulty unless a par- 
ticularly poor guess of the variables at Xm is employed. The 
magnetic diffusion shock is continuous and may be inte- 
grated through without pause. As magnetic diffusion is not 
important in this region the integration to the boundary 
(typically located at Xout = 10*-10^) is rapid and uncompli- 
cated. Integrating in the inwards direction is more problem- 
atic as the calculation is very sensitive both to the values of 
the variables at Xm and the Hall diffusion parameter. 



3.2 Inwards integration to Xc 

Integrating inwards is complicated by the presence of shock 
discontinuities and sonic points, which must be calculated 
explicitly. Close to the magnetic diffusion shock there may 
occur a subshock in which the supersonic (but slowing due 
to the sudden increase in fez) inflow is abruptly slowed to a 
subsonic rate. This only occurs when the Hall parameter is 
positive, and is likely caused by the changed magnetic brak- 
ing triggered by the azimuthal fleld growth during the shock. 
This subshock is a sharpening of the post-shock variation of 
the density and infall speed in the solutions of KK02. 

The magnetic diffusion subshock's existence is detected 
by performing a test integration: if the variables approach 
a sonic point (where = 1) then a shock must exist up- 
stream. Its position is found using the same method as that 
for the centrifugal shock described below, using the same 
jump conditions. Downstream of this subshock a sonic point 
occurs as the radial velocity becomes supersonic once more; 
a small manual step through the sonic point is performed 
and the variables are integrated to the centrifugal shock. 

The position of the centrifugal shock is found by per- 
forming a binary search over an appropriate interval. The 
upper and lower bounds on Xc are described by Xco ± 0.2a;co 
where Xco is estimated from KK02's equation 65, using 
(t^a - Sfin) in place of fjA- 



XcO 



rric exp 



23/2r 



(84) 



This is an overestimate of the true shock position (typically 
by around 20 per cent), however the routine is sufficiently 
robust that this does not present a problem. The variables 
are integrated to the estimated position of Xco, where the 
jump conditions derived below are applied, and then inte- 
grated towards the inner boundary. Unless the shock po- 
sition is known very precisely, the variables will approach 
their asymptotic values and then veer off course. 

This behaviour is most clearly seen in the surface den- 
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sity a, which increases rapidly downstream of the shock if 
Xco is an overestimate to the true shock position, and drops 
dramatically if Xco is an underestimate. The incorrect esti- 
mate is then assigned to be the new boundary of the shock 
position as appropriate, and a new estimate of Xco is chosen 
at the midpoint between the boundaries. As the position of 
the shock is more precisely known, the variables follow the 
asymptotic behaviour longer. 

The potential presence of any sonic points and sub- 
shocks does not interfere with the iterative routine for find- 
ing the shock position. When there exists a sonic point 
downstream of the shock: a bounces upwards at the sonic 
point if XcQ is too high, while if Xco is too low then the inte- 
gration fails at the sonic point. When the shock position is 
known to approximately 10"^", it is considered to be known 
to the precision of the full calculation. 

The variables are integrated inwards from the centrifu- 
gal shock and if rjH is large enough they are manually inte- 
grated through the sonic point that presages the existence 
of a centrifugal subshock. After the sonic point, the binary 
search is employed again to find the position of its associated 
shock front, as the behaviour of the variables downstream 
of the subshock is unchanged. In this instance the initial 
upper boundary of the search (a^up) is the sonic point, and 
the lower boundary is Xdown ~ O.lXup. If fjH is so large that 
multiple subshocks exist then this behaviour repeats, with 
the first subshock followed by an additional sonic point and 
then a second subshock; the increasing number of subshocks 
cause the calculation to become so unstable that it is unable 
to converge on the true similarity solution. 

Note that even when the shock positions are known 
"precisely", the routine is often unable to complete the in- 
tegration all the way to the inner boundary. This is avoided 
by integrating a simplified set of equations from a point 
far from the centrifugal shock to the inner boundary. These 
equations and the matching criteria are outlined in ^3.41 



integrated over the shock to give the relations 

aw — constant (85) 

and a (w'^ + = constant, (86) 

which apply when the magnetic forces are small. These equa- 
tions give the non-trivial jump conditions: 

wd = — (87) 
and ad = auW^, (88) 

where u and d denote the upstream and downstream sides of 
the shock; these are also used for any subshocks downstream 
of the principal shocks. KK02 also derived jump conditions 
for use when the magnetic forces are important; these are not 
required in this work. It is possible that as the positive Hall 
parameter becomes larger the magnetic diffusion subshock 
may require different jump conditions that take the twisting 
of the magnetic field lines due to the Hall effect into account. 

The magnetic field strength does not change in the cen- 
trifugal shock front as the magnetic pressure and tension 
terms are not large enough to change the field behaviour. In 
the post-shock region the field increases as the other vari- 
ables settle to the asymptotic behaviour in ij2.3l While the 
magnetic field is unchanged by the passage of the shock, 
its position and strength depend upon the Hall parameter. 
The increased magnetic diffusion can also cause the forma- 
tion of subshocks - rings of sharply-enhanced density in the 
post-shock region - as the outward-moving fiux causes the 
infalling gas to be slowed. 

As with the magnetic diffusion subshock, the centrifugal 
subshocks are preceded by a sonic point that must be cal- 
culated manually and the subshock position is found using 
the binary search above. Inwards of the centrifugal shocks 
the variables approach the asymptotic behaviour with some 
overshoots and corrections. 



3.3 Jump conditions 

The magnetic diffusion shock is smooth and continuous, 
and does not require any explicit calculation of shock con- 
ditions. The shock is the transition between the approxi- 
mations 6z,!ou) and bz,high (Equations [7S] and I79|) . and in 
the shock front b^, b^^s and h are changed, with fo^.s down- 
stream of the shock given by b^^s ~ —wxdSbz,iow During 
the shock the scale height is markedly compressed, suggest- 
ing that a breakdown of vertical hydrostatic equilibrium has 
occurred. In reality the enhanced magnetic squeezing during 
the shock front would be unable to reduce the disc thickness 
so dramatically over the fluid transit time through the shock, 
which implies that the thinness of the compressed region is 
a numerical artefact. Furthermore, as the magnetic pressure 
far exceeds the gas pressure any breakdown of isothermality 
would not greatly affect the collapse. 

The centrifugal shock is a true discontinuity in the 
fluid variables, which must be explicitly calculated. The 
shock conditions used are those for the "isothermal jump" in 
KK02; these are found by recognising that the derivatives of 
the surface density and radial velocity at the shock are large 
compared to the other terms in the conservation of mass and 
radial momentum equations (|43l and 144 |l . The equations are 



3.4 Innermost integration 

In the innermost regions of the collapse the derivatives 
da/dx and dbz/dx become very large in comparison to the 
other derivatives. These can correspondingly cause small nu- 
merical errors in the calculation of these derivatives and 
their integrals to build and trigger the appearance of a spon- 
taneous singularity where the variables diver ge dra matically 
from the asymptotic solution (discussed in lLilllQQ^ ). This be- 
haviour is avoided by replacing the full MHD equations with 
a simplifled set that can be integrated all the way to the in- 
ner boundary: the problematic derivatives are replaced by 
approximations derived from the expected asymptotic be- 
haviour, that is, 



and 



dbz 
dx 
da 
dx 



5 m. 
"i 

3 
'2 



3/4 



-9/4 



25 



-5/2 



2V(2<5//)^ + 1 



(89) 



(90) 



These are then substituted into the other equations so that 
a, bz, b^^s and h are found by solving Equations 1441 [46l 1471 
and simultaneously, while the remaining equations are 
integrated to the boundary. 

Switching between the full and simplifled models occurs 
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Table 1. Boundary conditions and parameters. 



parameter 




value 


ambipolar diffusion 


VA 


1 nn 
i.UU 


cap on b^^s 


S 


1 on 


magnetic braking 


a 


0.08 


boundary condition 


value 


mass-to-fiux ratio 


MO 


2.90 


rotational velocity 


vo 


0.73 


radial velocity 


UQ 


-1.00 


core accretion rate 


A 


3.00 



after the minimum in the surface density as the variables fol- 
low the asymptotic solution, but before they diverge from 
the expected behaviour. Typically this occurs when the old 
and new values of a match to within 0.01/a; (~ 0.1 per cent 
of the original value) and da/dx calculated using both meth- 
ods matches less well to around 200 /x (~ 7 per cent). When 
the precision of the match is as high as possible the change 
between the two equation sets is smooth; the required pre- 
cision of the match is raised when the switch is apparent as 
this clearly influences the accuracy of the calculations. 

The simplified set of equations is not subject to the 
same numerical instabilities as the full set unless the guess at 
the matching point is particularly poor. Should such a guess 
be adopted then the routine is unable to match onto the 
simplified model, and the failure of the inwards integration 
is used to refine the values at the matching point. 

The integration is continued to the inner boundary (typ- 
ically at Xin = 10~*) where the variables are compared to 
the inner conditions given by Equations 1641 and 1711 The dif- 
ferences between the expected and integrated variables are 
passed back to the shooting routine, which modifies the val- 
ues at Xm and begins the next iteration. The similarity so- 
lution is considered to have converged when the integrated 
variables match the boundary conditions to at least 0.002 
per cent. Due to the ease of convergence in the outwards di- 
rection, the outer variables match to a much higher degree. 



Ambipolar Diffusion Coiiapse 
10"^ 10-2 10"'' 10° lO"! 10^ 




Figure 1. Calculated similarity solution for collapse with only 
ambipolar diffusion, duplicating KK02's fig. 7. The enclosed mass, 
accretion rate, surface density and scale height are displayed in 
the upper panel; the magnetic field components and enclosed flux 
in the centre panel; and the gas velocities and angular momentum 
in the lower panel. The parameters and boundary conditions are 
as in Table [T] the central mass is rric = 4.67; and the magnetic 
diffusion and centrifugal shocks are located at Xd = 0.406 and 
Xc = 1.32 X 10~^ respectively. 



4 RESULTS 

Superficially, the similarity solutions with Hall diffusion are 
similar to the fiducial ambipolar diffusion solution of KK02. 
We match their boundary conditions and parameters, listed 
in Table [T] holding the ambipolar diffusion parameter con- 
stant at rjA = 1-0 so that the only changes between the sim- 
ilarity solutions presented in their work and ours are those 
wrought by the addition of Hall diffusion. 



4.1 Ambipolar diffusion solution 

As a test of the computational code we calculate the fiducial 
ambipolar diffusion solution of KK02 in Fig. [T] and follow 
their discussion of the solution here. The outer regions of all 
the solutions match IMHD similarity solutions, as the mass- 
to-flux and mass-to-angular momentum ratios are constant 
while the material falls in at supersonic speeds. The radial 
velocity and scale height are dominated by the self-gravity of 



the disc, which causes the material to fall towards the mid- 
plane before it is accelerated towards the central mass. The 
magnetic fleld gradually builds up as the matter falls inward, 
becoming important to the dynamics at around x ~ 2 where 
magnetic braking starts to affect the angular momentum 
transport and the constant ratio of the mass-to-angular mo- 
mentum breaks down. The azimuthal fleld attains its capped 
value, and becomes important to the angular momentum 
transport. 

The mass and angular momentum in this region tend 
towards plateau values: the mass to that in Equation 1761 
and the angular momentum to 



(91) 



which depends upon the boundary conditions. This plateau 
forms as the gravity of the central point mass starts to dom- 
inate and the gas starts to fall inward rapidly. The surface 
density and magnetic field build up and ambipolar diffusion 
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Figure 2. Schematic of tiie poloidal magnetic field in tlie mag- 
netic diffusion sfiock. The disc (dotted Unes) is compressed as 
the vertical field becomes large, causing the field lines at the sur- 
face to straighten from being largely radial upstream to having 
roughly equal poloidal components downstream. (Not to scale.) 

becomes important, causing flux-freezing to break down at 
the magnetic diffusion shock (at Xd = 0.406). 

This shock takes the form of a sudden increase in the 
vertical magnetic field as the field lines diffusing against 
the flow in the downstream ambipolar diffusion-dominated 
regime meet those coming inward under IMHD collapse. The 
field geometry is illustrated in Fig. [2l upstream the surface 
magnetic field is dominated by the radial component, which 
is an order of magnitude larger than the vertical and az- 
imuthal components. During the shock (and the downstream 
transition region) the field lines straighten until the poloidal 
components are approximately equal at the disc surface. The 
increase in field strength slows the gas, and magnetic squeez- 
ing comes to dominate the vertical compression. 

Immediately interior to the magnetic diffusion shock 
the poloidal field components scale with ~ x^^ as the 
surface density and thickness of the disc increase and the 
gas is slowed by the radial magnetic pressure. Rotation 
is not dynamically important and the shock has a similar 
structure to that in the nonrotating similarity solutions of 
IContopoulos et al] (|l998f ). The gravity of the central mass 
becomes the dominant radial force on the gas at the end of 
the post-shock region (the peak in the surface density) and 
the matter is accelerated inwards until it is in near-free fall 
collapse. 

As in slo wly-rotating sim ilarity solutions that are non- 
magnetic (e.g. lSaigo fc Hana wa 1998) or include IMHD (e.g. 
KK02 §3.2) rotation remains dynamically unimportant un- 
til the vicinity of the centrifugal shock. Between the two 
shocks the gravity of the central mass dominates the radial 
infall of the gas, which is slowed only a little by ambipolar 
diffusion of the magnetic flux and the magnetic pressure. 
The magnetic braking increases the azimuthal field com- 
ponent until it is again capped; the activation of the cap 
causes the change in the behaviour of the scale height at 
X ~ 4.5 X 10~^ as 6^,3 then contributes more to the verti- 
cal squeezing forces. The enclosed mass and accretion rate 
flatten and remain near-constant throughout the remainder 
of the collapse. The angular momentum starts to plateau 
again and the centrifugal force becomes large and equal to 
the gravitational force, triggering the centrifugal shock at 
Xc = 1.32 X 10"^. 

The centrifugal shock slows the gas so that the infall is 
subsonic and the surface density increases by more than an 



order of magnitude. The shock is followed by a thin layer in 
which the azimuthal and vertical magnetic field components 
increase rapidly. This causes the angular momentum to drop 
to its asymptotic behaviour, as the surface density and the 
other variables adjust with a few overshoots towards their 
expected rotationally-supported disc behaviour. The transi- 
tion between the full model and the simplified set of equa- 
tions occurs at a; ~ 8.6 x 10~^, after the variables have joined 
onto the asymptotic disc described by Equations I64I[7^ 

The Keplerian disc itself is small, with a mass ~ 5 per 
cent that of the central point mass. The nondimensional 
mass at the origin is rric = 4.67, corresponding to a moderate 
accretion rate of Mc ~ 7.6 x 10~® M© yr~^ (so that at a time 
t = 10* yr, the central mass is Mc = 0.76 Mq). The surface 
density of the disc depends upon the ambipolar diffusion and 
azimuthal field cap parameters, as does the infall velocity, 
which is subsonic and very low. The disc is extremely thin, 
and the vertical squeezing is dominated by the tidal and self- 
gravitational forces. The magnetic flux in the disc scales as 
i'^^'*, so that — >■ as a; — >■ 0; clearly the amount of flux 
present in the protostar depends upon more detailed flux 
transport and destruction mechanisms t han are included in 
this model, such as Ohmic diffusion (e.g. Li & McKee"l996|) 
and reconnection (e.g. iGalli fc Shuil993b. : .Lazarian.i2005i ). 

4.2 Hall diffusion solutions 

The flrst of the similarity solutions with Hall diffusion is that 
presented in Fig. [3]for the self-similar collapse of a molecular 
cloud with fjH ~ —0.2 and other parameters matching those 
in Fig. [T] The negative Hall parameter solutions have more 
radial diffusion of the magnetic fleld against the neutral fluid 
and charged grains, so that the magnetic pressure builds up 
earlier in the collapse process, triggering the formation of 
the magnetic diffusion shock. The negative Hall parameter 
also increases the initial rate of magnetic braking so that 
&0,s attains its capped value earlier in the collapse, and the 
magnetic braking is then determined by the strength of the 
vertical field component. 

As in the no-Hall solution, at the outer edge of the col- 
lapse in Fig. [3] the matter is falling in supersonically under 
IMHD. As the surface density builds up the fleld does too, 
causing the magnetic pressure and tension terms to become 
important, while the magnetic braking transports angular 
momentum from the infalling gas to the external envelope. 
The angular momentum and enclosed mass start to plateau 
as the dominant force on the radial velocity switches from 
the self-gravity of the disc to the gravity of the central mass, 
which in turn causes the accretion rate to taper off. The 
formation of the magnetic diffusion shock at Xd = 0.461 
(increased from the non-Hall solution) is caused by the de- 
coupling of the field from the neutral particles. 

The magnetic diffusion shock in this solution is weaker 
than in Fig. [T] as most of the neutral particles and grains 
have already decoupled from the magnetic field, so that the 
vertical field component increases by only 4.5 times (c/. 6.2 
times in Fig.[TJ. The disc is less vertically compressed by the 
fleld, producing a thinner shock, and b^^s does not grow as 
rapidly. Within the shock the fleld is further decoupled from 
the neutrals and grains, allowing Hall and ambipolar diffu- 
sion to become more important downstream of the shock 
and throughout the remainder of the solution. The field lines 
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Figure 3. Hall collapse with f]ii = —0.2. The variables are as in 
Fig-El the parameters and boundary conditions are given in Table 
[T] The central mass is rac = 4.23; the magnetic diffusion and 
centrifugal shocks are located at = 0.461 and Xc = 3.78 X 10"^ 
respectively; these are increased from the non-Hall positions, and 
the post-shock regions are smoothed by Hall diffusion. 
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Figure 4. Collapse with larger Hall parameter f]ii = —0.5. The 
boundary conditions and parameters match those in Fig. |3] (Table 
[TJ. The central mass is reduced to rric = 3.77 and the magnetic 
diffusion and centrifugal shocks moved outwards to x^j = 0.557 
and Xc = 8.31 X 10~^ as the increased magnetic diffusion smooths 
the post-shock regions and increases the size of the Keplerian disc. 



straighten as in Fig. [J] although the radial field component 
is still dominant, the vertical component increases in the 
shock until it is just smaller than br,s- 

Downstream of the magnetic diffusion shock the surface 
density gradually increases as the infall velocity is slowed 
by the larger magnetic support. This post-shock region is 
smoother than that without Hall diffusion, presenting a 
gentler transition to the free fall collapse that occurs out- 
side of the rotationally-supported disc. The vertical field 
scales as in this region as the increased radial diffusion 
means that there are fewer field lines in total moving against 
the fiow of the neutral particles. The magnetic braking de- 
creases the angular momentum efficiently until b^^s attains 
its capped value and j begins to plateau once more. 

The centrifugal force builds up and triggers the cen- 
trifugal shock at Xc = 3.78 x 10"^ (c/. Xc = 1.32 x 10"^ 
in Fig. [l|. The shock is a discontinuity in the surface den- 
sity and radial velocity, which is again less strong than in 
the solution without Hall diffusion, and inwards of this the 
vertical and azimuthal field components increase steeply as 
the field reacts to the shock (although b^^s remains capped 



at —5bz). Downstream the variables tend (with overshoots) 
towards their asymptotic values. 

The inner disc is in Keplerian rotation satisfying Equa- 
tions [64HZS] The central mass is nic = 4.42, decreased from 
the non-Hall solution, and corresponds to an accretion rate 
of Mc = 7.21 X 10"*^ Mo yr"\ The scalings of the other 
variables with respect to x are the same as in Fig. [1] how- 
ever the surface density is increased by the larger magnetic 
diffusion parameter / = 1.72 (c/. fifjH = 0) = 1.3); and 
the increased radial magnetic diffusion causes the strength 
of the magnetic field to be decreased from that in the am- 
bipolar diffusion-only solution. This in turn means that less 
matter can lose its angular momentum and fall onto the cen- 
tral mass, so that the gas is at a higher surface density in 
this larger Keplerian disc. 

The next similarity solution, presented in Fig. |4l shows 
the calculation with j)h ~ —0.5 on the same scale and with 
the same parameters as Fig. |3l The total radial magnetic 
diffusion is further increased so that many of the neutral 
particles have decoupled from the field before the magnetic 
diffusion shock at Xd ~ 0.557; this causes the intensity of 
the shock to drop further so that the vertical field strength 
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Figure 5. Gravitational collapse with positive Hall parameter 
fiH = +0.2. The boundary conditions and parameters match Fig. 
[3] and Table [T] the central mass is rric = 4.63. The positive Hall 
term causes the formation of subshocks: the magnetic diffusion 
shocks are located at a;^ = 0.365 and x^2 = 0.260; the centrifugal 
shocks at Xc = 6.05 X lO"'^ and Xc2 = 5.21 X 10^'^. 



is only 4 times larger than its original value. There is less 
of a magnetic wall at this point as less flux remains to be 
decoupled from the neutrals within the shock itself. 

As in the previous solution with negative Hall param- 
eter the post-magnetic diffusion shock region is smoothed, 
with even less change in the surface density and radial ve- 
locity. The gas is slowed by the magnetic diffusion shock, 
but the gravity of the central mass quickly overcomes this 
and pulls the fluid inwards. The radial velocity downstream 
of the post-shock region increases as the fluid nears the pro- 
tostar, however it remains below the free fall velocity at all 
times. The mass and angular momentum both plateau in 
this region before the increasing centrifugal force triggers 
the centrifugal shock. 

The centrifugal shock occurs much earlier in this sim- 
ilarity solution at Xc = 8.31 x 10~^. This is brought about 
by the decreased values of bz and fc^.^ in the free fall region, 
which reduce the amount of magnetic braking that takes 
place and cause the centrifugal force to become important 
earlier. Inwards of this shock is a much wider region of ad- 
justment as the variables join the inner disc solution. 

The Keplerian disc is substantially larger than that in 




Figure 6. The magnetic diffusion shock and subshock Xd2, for 
the collapse with fj[f = -1-0.2. As seen in Fig. \E\ j and are 
continuous, and —u mirrors the behaviour of crx. 



the previous solution, containing ~ 38 per cent of the mass of 
the central protostar. The surface density has also increased 
as the magnetic diffusion parameter is / = 2.31, while the 
lowered central mass rric = 3.77 corresponds to a central ac- 
cretion rate of Mc = 6.15 x 10~® Mq yr~^. Again, the larger 
disc corresponds to a lower accretion rate, as the reduced 
magnetic braking prevents the fluid from losing rotational 
support and falling in. 

The flnal solution presented here is that with fjH ~ +0.2 
in Fig. [51 which is the most dynamically different from those 
of KK02. Although the initial conditions and parameters 
match those in Fig. O the change in the sign of the Hall 
parameter, which corresponds to a reversal of the orientation 
of the magnetic field with respect to the direction of rotation, 
introduces many changes to the collapse dynamics. 

These begin at the magnetic diffusion shock, which has 
moved inwards from f]H = solution to Xd ~ 0.366. This 
shock is of increased intensity due to the reduced radial mag- 
netic diffusion upstream, which causes a larger increase in 
in the shock. The magnetic braking downstream is increased 
by the presence of a stronger field and the sign of the Hall 
term in Equation 1521 The disc is more sharply compressed 
as the field lines straighten at the shock front, and the fluid 
is so slowed by the increase in the magnetic pressure that a 
second shock front forms at Xd2 = 0.260. 

In the magnetic diffusion subshock, shown at higher res- 
olution in Fig. the fluid is slowed until the radial veloc- 
ity is low and subsonic. The surface density increases under 
the jump conditions from i]3.3l this ring of matter contains 
approximately 18 per cent of the protostellar mass. The az- 
imuthal fleld component and the scale height also increase, 
while dbz / dx decreases steeply. The infall region downstream 
of the subshock is wider in logarithmic similarity space than 
in the previous solutions, with the increased magnetic brak- 
ing reducing the angular momentum more quickly as the 
radial velocity increases. The surface density drops as the 
fluid falls in and magnetic squeezing dominates the vertical 
compression until the gravity of the central mass takes over 
near the centrifugal shock. 

The centrifugal shock occurs at Xc = 6.05 x 10""^, half 
that of the fju = solution, dramatically decreasing the size 



© 2012 RAS, MNRAS 000.[Tlf24l 



14 C. Braiding and M. Wardle 



10° 



1 1 1 1 


- 


! Zs'-'' / ^ 




lOOh/x 

: xb^ 




Oj: 

I.I. 



4x10-3 5x10-3 6x10-3 7x10-3 
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Figure 8. The centrifugal shock radius (in au at IC* years) 
against the ratio of the Hall to ambipolar diffusion parameters. 



of the rotationally-supported disc. This change is brought 
about by the increased magnetic braking caused by the pos- 
itive Hall parameter, which reduces the angular momentum 
so that the centrifugal force is not dynamically important 
until the gas is very close to the protostar. Downstream of 
the shock the slowed fluid accelerates inwards as the mag- 
netic field increase forces additional magnetic braking and 
a drop in the centrifugal force. The surface density drops as 
the infall velocity becomes supersonic, and the centrifugal 
force becomes important once more, triggering the subshock 
at Xc2 = 5.21 X 10""^ which is shown in Fig.[71 The matter is 
slowed in the subshock until it is again subsonic. Although 
the surface density increases in the subshock, the disc does 
not beco me gravitation ally-unstable as the Toomre Q pa- 
rameter (|Toomrelll964h remains above 8 in this region. 

Downstream of the subshock the variables again settle 
with overshoots to the Keplerian disc behaviour. The central 
mass is rric — 4.63, corresponding to a protostellar accretion 
rate of Mc = 7.53 x lO"*^ yr-\ while the disc con- 
tains only 1.6 per cent the protostellar mass. The magnetic 
diffusion parameter is / = 0.945, reflecting the decrease in 
surface density as the magnetic field strength increased. 

The magnetic diffusion and centrifugal subshocks both 
occur only in the calculations where fjH is positive, as the 
larger magnetic pressures and braking caused by the increase 
in the infalling magnetic field force the gas to rapidly change 
in radial velocity and density. The number of centrifugal sub- 
shocks increases with fjH — three subshocks have been ob- 
served in one similarity solution that was not properly con- 
verged (due to the numerical instability of the subshocks) at 
the time of publication — while only one magnetic diffusion 
subshock has been observed. 

The diffusion parameter of the disc, /, cannot be less 
than or equal to zero in our collapse solutions, which do not 
contain counter-rotation. This restriction limits the range 
of positive r)H that can be explored: as fjH increases, the 
size and surface density of the rotationally-supported disc 
decrease, and the rings of gas formed by the subshocks 
are more likely to be gravitationally unstable. Similar con- 
straints were found to limit the laun ching of disc winds in 
the analysis of lSalmeron et all l|201ll '). who showed that disc 
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Figure 9. The protostellar accretion rate (in 10"® M0/yr) 
against the ratio of the Hall to ambipolar diffusion parameters. 



wind solutions only exist for particular combinations of the 
field polarity and the ratio of the Hall to ambipolar diffusion 
parameters. 

The positions of the shocks change with the Hall param- 
eter as demonstrated in Fig. [S] which plots the dimensional 
centrifugal shock position for similarity solutions with pa- 
rame ters as in Table [T] at 10^ years (illustrated in iBraidind 
l201lf ) against the ratio of the Hall to ambipolar diffusion 
parameters. While both directions of Hall drift contribute 
to the size of the disc, the radial drift of (which increases 
when rjH < 0) has a greater effect on the radius of the cen- 
trifugal shock than the azimuthal Hall drift. 

There also exists a correlation between the accretion 
rate onto the central protostar and t he radial magnetic 
field diffusion, shown in Fig. [9] (see also IContopoulos et al] 
Il998l ). The accretion depends upon the disc radius, with 
larger discs correspondin g to lower accretion rates and vice 
versa (| Allen et al.ll2003bl ). as a more negative Hall param- 
eter causes increased drag on the neutrals and reduces the 
radial velocity of the fiuid. The disc radius also depends on 
the initial rotational velocity, as the centrifugal force is more 
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important and a larger disc forms when the initial angular 
momentum of the core is large (KK02). 

The protostellar accretion rate appears to turn over 
at around Mc ~ 7.6 x 10~® Mq yr~^ as the ratio of the 
Hall to ambipolar diffusion parameters becomes positive and 
greater than 0.05. This is due to the formation of subshocks 
in the solutions with fjH > 0: as the density is enhanced 
accretion through the disc drops. As the Hall parameter in- 
creases further subshocks are introduced, the diffusion pa- 
rameter tends towards zero and the rings formed by the 
subshocks become more unstable. 



5 DISCUSSION 

The similarity solutions clearly show that Hall diffusion 
changes the structure and dynamics of the collapse of molec- 
ular cloud cores into protostars and protostellar discs. The 
rotationally-supported disc size, and the accretion rate onto 
the protostar are determined by the ratio of the Hall and am- 
bipolar diffusivities, which influences the magnetic braking 
affecting the rotation of the collapsing core. It is also clear 
that Hall diffusion can inhibit disc formation by enhancing 
the magnetic braking, or by counteracting ambipolar diffu- 
sion to the point that the field starts to infall faster than 
the fluid, increasing the magnetic pressure and tension. 

The dependence of the similarity solutions on the orien- 
tation of the magnetic field and the sign of the Hall diffusion 
parameter f]H (more specifically upon the sign of fini^ ■ fi)) 
gives rise to two different patterns of collapse behaviours. 
The similarity solutions with f}H = Q and ±0.2 are con- 
verted to dimensional form and plotted against the radius r 
(at t — 10* years) in Fig. 1101 with the surface density in the 
upper panel and the vertical magnetic field strength plotted 
as in the lower. The solutions have the same boundary 
conditions and parameters (Table[T|, and the surface density 
and vertical field strength (at r = 1 an), the central mass, 
and the mass and size of the inner disc (all a± t — 10'* years) 
are listed in Table [21 The outer regions where IMHD holds 
are near-identical and it is only near the magnetic diffusion 
shock at r ~ 100 au that the changes brought on by Hall 
diffusion become apparent. 

The dotted lines in Fig. 1101 corresponding to the neg- 
ative Hall solution in Fig. [S] show the formation of a large 
rotationally-supported disc that has radius Rc ~ 15 au at 
(t = 10* years) and the highest inner Keplerian disc surface 
density of all the solutions. The dashed lines are the f]H = Q 
solution from Fig. [T] which possesses a disc radius half that 
the negative Hall solution. The surface density of the Kep- 
lerian disc has decreased by a constant factor from that in 
the negative Hall solution. Finally, the solid curves charac- 
terise the similarity solution with positive Hall parameter 
f]H = 0.2 (Fig. [5]) which has a Keplerian disc that is almost 
an order of magnitude smaller than that in the negative Hall 
case. This disc is bounded by a thin ring of enhanced den- 
sity that rapidly drops off as the magnetic field peaks; the 
material is then shocked again and comes to match onto the 
inner solution. The density is much lower than in previous 
solutions, and the disc grows at a slower rate. 

The similarity solutions span many orders of magnitude 
in both radius and density, and the inclusion of a Hall pa- 
rameter that is 20 per cent that of the ambipolar diffusion 



10"' 10° 10^ 10^ 10^ lO'* 



10° 









flH = -0-2 








— -?)H = ■ 








fi„ = +0.2 










; 1 1- 


1 




1 








— -nH = o 








— •nH = +0-2 ; 


- 

















-21 Hq- 

10"'' 10° 10' 10^ 10^ lO'* 



r @t= 10000 yrs (au) 

Figure 10. The surface density E and the vertical magnetic field 
component Bz plotted against radius at t = 10* yr for the solu- 
tions with fiA = 1.0 and fiu = —0.2 (dotted line), (dashed line) 
and -1-0.2 (solid line). These solutions were plotted individually 
in nondimensional form in Figs. |3] [Hand [5] respectively. 

parameter has a large effect on the behaviour of the mag- 
netic field. In the intermediate region between the magnetic 
diffusion shock (at Rd ~ 100 au in Fig. llOp and the centrifu- 
gal shock [Rc ~ 1 — 10 au), the azimuthal field tension causes 
the Hall drift to enhance the radial diffusion of the field lines 
when fjH is negative. The magnetic diffusion shock occurs 
earlier in the collapse and is less dynamic than in the other 
solutions, as much of the field has already been decoupled 
from the fluid. 

However, when f}H is positive then Hall diffusion acts to 
reduce the net radial diffusion, resulting in magnetic walls 
and subshocks that disrupt the flow. The magnetic field car- 
ried inwards increases, and the magnetic pressure and ten- 
sion terms remain important throughout the collapse. Any 
twist in the field lines causes an increase in the magnetic 
pressure gradient, so that the net amount of radial diffusion 
drops off as the magnetic braking slows the rotation. 

There is a similar duality to the azimuthal field drift. 
Again looking at the region between the two shocks, when 
fjH is negative Hall drift occurs in the azimuthal direction, 
twisting up the field lines in the pseudodisc and creating a 
leading torque on the neutral rotation. The reduced value of 
Bz causes the azimuthal field component to reach its capped 
value -B</,,s = —5Bz sooner, and the magnetic braking, which 
depends upon BzB^^s, is also reduced. Because of this less 
angular momentum is removed from the pseudodisc, caus- 
ing the centrifugal force to become dynamically important 
earlier and a larger rotationally-supported disc to form. 

In the other orientation when fjn is positive. Hall and 
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Table 2. The surface density and vertical field component in the Keplerian disc at r = 1 au, the protostcUar 
mass and the size and mass of the Keplerian disc at t = 10'* years for the solutions depicted in Fig. 1101 



Vh 


S (g cm ^) 


(G) 


M, (Mq) 


Mmsc (Mq) 


Rc (au) 


-0.2 


1920 


0.289 


7.21 X 10-2 


9.99 X 10-3 


15.10 





1250 


0.299 


7.62 X 10-2 


3.75 X 10-3 


5.31 


0.2 


620 


0.304 


7.54 X 10-2 


1.24 X 10-3 


2.43 



ambipolar diffusion act together to untwist the field lines in 
the pseudodisc. In these similarity solutions is larger, and 
so while it takes longer for B^.s to achieve its capped value 
there is more magnetic braking and the angular momentum 
is further reduced. A smaller Keplerian disc forms due to the 
reduced centrifugal force, and both shocks have subshocks 
where the magnetic forces alter the radial velocity of the 
fluid. Downstream of the magnetic diffusion shock the radial 
magnetic pressure gradient slows the fluid in the radial and 
azimuthal directions, while downstream of the centrifugal 
shock the gas is accelerated inward as the increase in Bz 
causes a burst of magnetic braking that disrupts the disc 
and causes the formation of a subshock. 

While the angular momentum behaviour between the 
shocks is changed by the inclusion of Hall diffusion, the 
cap on B^^s acts to ensure that the angular momentum in 
the inner disc is that expected for a Keplerian disc. The 
cap, while physically motivated, replaces unspecifled disc 
physics such as reconnection, a disc wind, or turbulence, 
which would act to prevent the azimuthal fleld component 
from greatly exceeding the vertical component; and the mag- 
nitude at which it ought to act to limit B^^s is uncertain. 
It is also unclear if such limiting of the azimuthal fleld com- 
ponent happens in real collapsing cores, as numerical sim- 
ulations do demonstra te tightly-wound magnetic flelds (e.g. 
iMachida etHI l2008bl ). The azimuthal field cap limits the 
similarity solution set explored to those in which discs form, 
however despite this Hall diffusion has been shown to re- 
strict disc formation in collapse without counter-rotation if 
the Hall diffusion is too strong in comparison to ambipolar 
diffusion and t^h has the "wrong" sign. If the core was ini- 
tially rotating in the opposite direction, the opposite sign of 
fjH would be problematic to disc formation. 

All of the solutions in Fig. [10] and iBraidind (|201ll ) 
form protostars of around 0.7 solar masses with protostellar 
discs of radius Rc ~10-150 au and mass Md ~ lO-^-lO-* 
Mq in t = 10^ years; these are the same order of mag- 
nitude expected from observations of Class I YSOs (e.g. 
|j0rgensen et al.l |2007| ) . The surface density of the disc is 
quite sensitive to the Hall diffusion parameter, and scales 
as E oc r-3''2 and E cx t^^^ in the inner Keplerian disc, with 
values E(r = 1 au, t = 10* years) ~ 10^ g cm-^ (see Table 
[2|. These values of the surface density are consistent with 
what is thought to have occurred in the solar nebula (e.g. 
the mini mum mass solar nebu la has S — 1700 g cm-^ at 
r = 1 au; IWeidenschillin3ll977f) . 

It has recently been argued that there exists a hand- 
edness to observations of transverse gradients in the Fara- 
day rotation measure across the base of jets associated with 
active galactic nuclei (AGN: [Contopoulos et al]|2009l '). The 
majority of sources in which it was possible to determine 
the transverse gradients were found to have clockwise gradi- 
ents, implying that the outflow has a helical magnetic fleld 



with a preferred magnetic polarity. One explanation of this 
behaviour is that the Hall effect is important in the inner ac- 
cretion disc, acting to form a jet when the fleld has a positive 
magnetic polarity, a nd to suppre ss jet formation when the 
polarity is negative (|Konigjl2010l ). This explanation flts the 
limited available data well, although it must be conflrmed by 
future observations at higher resolutions and sensitivities. 

Similarly, it may be possible to show the importance 
of the sign of the Hall parameter observationally by mea- 
suring the polarisation of the magnetic fleld with respect to 
the axis of rotation in Zeeman observations of newly-forming 
stars and their discs. Should larger discs and lower accretion 
rates be correlated with a particular fleld orientation then 
the Hall effect will have been shown to affect the collapse 
process. ALMA (among other next-generation instruments) 
shall be capable of imaging nearby dense prestellar cores and 
their envelopes in both dust and molecular line emission, and 
could also be used to observe polarised dust emission and 
map the magnetic fleld in cores. Such observations could 
be sensitive enough to observe if there is any difference in 
the fleld alignment between protostellar discs and their en- 
velopes, and whether there is any correlation between disc 
size and the direction of rotation in the disc. 

None of the solutions calculated the effects of very weak 
magnetic diffusion or strong magnetic braking on the core. 
In numerical simulations with such conditions the magnetic 
braking removes all of the angular momentum from the fluid, 
preventing the formation of a rotationally-supported disc 
this behaviour has been dubbed the "magnetic brakin; 



catastrophe" (see e.g. Allen et al.l 2003bl: 
l2009l : iHennebelle fc Ciardi|[2009l : iLi et al.l 



Mellon fc Umo. 



20111 '). Disc forma- 



tion was assured in this work by the cap placed upon B^^s, 
which limits the twisting of the fleld lines, however KK02 
were able to demonstrate the magnetic braking catastrophe 
by adopting a large 5. Hall diffusion i s capable of inducin g 
spin in an initially- nonrotating fluid (| Wardle fc NellT999l '). 
and could resolve the magnetic braking catastrophe in solu- 
tions with strong braking by spinning up the core in the op- 
posite direction once magnetic braking has removed the ini- 
tial angular momentum. Hall-induced s pin-up of an initially- 
nonro tating core was demonstrated by iKrasnopolskv et al.l 
({ 2011 ^, although the formation of a Keplerian disc required 
a larger Hall parameter than that expected from the micro- 
physics, possibly because the 10 au sink particle prevented 
the outwards growth of a disc from the geometric centre of 
the collapse. 

Further work must be done to study the role of the 
Hall effect on the magnetic fleld diffusion in star formation, 
particularly using the semianalytic model constructed in this 
work. This could include adopting more realistic values and 
scalings of the ambipolar and Hall diffusion parameters, in 
order to resolve the magnetic braking catastrophe; including 
a disc wind in place of the cap on B^,a; and exploring those 
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regions of parameter space in which Hall diffusion is the 
dominant form of field transport. 

As explained in i]2.1l the Hall and ambipolar diffusion 
coefficients in our calculations scale with the nondimensional 
variables in the same manner. However, the diffusion coef- 
ficients could scale as any function of x and the self-similar 
variables in order to mimic the behaviour expected from 
ionisation equilibrium calculations. For example, the Hall 
parameter could scale as 



(92) 



which is expected if the only particles in the collapsing core 
are neutrals and ions, without grains. Such a formulation 
would cause Hall diffusion to become important earlier in the 
collapse, and the Keplerian disc would need to be described 
by a new set of asymptotic inner boundary conditions. 

The inner Keplerian disc in our solutions satisfies the 
criterion Br.s/Bz > l/\/3, which is the launching condition 
for a cold, centrifugally-driven wind, and the radial scal- 
ing of the magnetic field components is identi cal to that of 
the ra dially-self-similar wind solution of .Bla n dford fc Pavnd 
||1982D . Such a disc wind (described in appendix C of KK02) 
would be the dominant mechanism for the vertical transfer of 
angular momentum from the disc to the envelope, and must 
be included in future self-similar collapse simulations in or- 
der to explore the influence a wind may have on the angular 
momentum transport and the magnetic braking catastro- 
phe, and to improve the accuracy of the semianalytic mod- 
els, as disc winds and j ets occur in numerical simulations 
of collapsing cores (e.g. IXomisakal I2OO2I : iMellon fc lH l2009l : 
ICiardi fc Hennebellell2010l l. some of which display the mag- 
netic braking catastrophe. 

Further explorations of parameter space are required to 
more fully understand the influence of the Hall effect on star 
and Keplerian disc formation. In particular we have yet to 
fully separate the influence of Hall diffusion from that of 
increased ambipolar diffusion in the core, to find similarity 
solutions in which Hall diffusion is the dominant form of 
flux transport, or solutions in which there is no ambipolar 
diffusion at all. As in iKrasnopolskv et al.l (|201lh . solutions 
with only Hall diffusion see the fleld diffusion depend upon 
the J X B terms, so the radial diffusion is controlled by the 
azimuthal fleld component (which may be capped), and the 
azimuthal diffusion depends upon the radial component. Ad- 
ditional work should also be done to confirm that larger discs 
form in Hall similarity solutions where the core is initially 
rapidly-rotating. 



6 CONCLUSIONS 

This paper described a semianalytic self-similar model of the 
gravitational collapse of rotating magnetic molecular cloud 
cores with both Hall and ambipolar diffusion, presenting 
similarity solutions that showed that the Hall effect has a 
profound influence on the dynamics of collapse. The solu- 
tions satisfled the vertically-averaged self-similar equations 
for MHD collapse under the assumptions of axisymmetry 
and isothermality, matching onto self-similar power law re- 
lations describing an isothermal core at the moment of point 



mass formation on the outer boundary and a Keplerian disc 
on the inner boundary. 

The inner solution describes a Keplerian disc in which 
accretion through the disc depends upon the magnetic diffu- 
sion; with an appropriate value of the nondimensional Hall 
diffusion parameter f]H sl stable rotationally-supported disc 
forms in which the surface density E scales as r~^^'^ and 
vertical field strength oc r~^^'^. These are the scalings ex- 
pected from other simulations of protostellar discs to which 
the solutions calculated in this work compare favourably. No 
disc may form in solutions without counter-rotation when 
the Hall parameter is large (in comparison to the ambipolar 
diffusion parameter) and has the wrong sign (which indi- 
cates the orientation of the magnetic field with respect to 
the axis of rotation) , as the diffusion in these solutions is too 
strong and causes disruptive torques that form subshocks in 
the similarity solutions. This behaviour occurs because the 
response of the fluid to Hall diffusion is not invariant under 
a global reversal of the magnetic fleld. 

The size of the rotationally-supported disc in the full 
similarity solutions was shown to vary with the amount 
of Hall and ambipolar diffusion affecting the pseudodisc 
through their effect on the magnetic braking in the fluid. By 
creating an additional torque on the disc. Hall diffusion can 
either increase or decrease the angular momentum and ro- 
tational support in the infalling fluid, leading to an order of 
magnitude change in the Keplerian disc radius between the 
similarity solutions at the extremes of —0.5 ^ ffH/fjA ^0.2 
(where the ambipolar diffusion parameter, t^a = 1) . A small 
amount of Hall diffusion was shown have a large effect on 
the solution because the dynamic range of collapse is itself 
many orders of magnitude in space and time. Hall diffu- 
sion causes there to be a preferred handedness to the fleld 
alignment and the direction of rotation in forming a large 
Keplerian disc that could be observed using next-generation 
instruments such as ALMA. 

The accretion rate onto the central point mass is sim- 
ilarly influenced by Hall diffusion. This is a smaller effect 
than that on the disc radius, as between fjn = ±0.1?7a (again 
with f}A = 1) the accretion rate onto the protostar only 
changes by 6 per cent, or 0.2 x 10~® M© yr~^. There exists 
a clear trend in which the protostellar accretion rate drops 
off with increasingly negative Hall parameter despite the 
constant accretion onto the core, as the reduced magnetic 
braking in these solutions causes a larger Keplerian disc to 
form, and accretion through this disc onto the protostar is 
slow. 

The magnetic braking catastrophe could be resolved by 
the inclusion of Hall diffusion in numerical solutions, as with 
one sign of f}H the Hall effect acts to reduce the total amount 
of braking affecting the core, preventing it from removing 
too much angular momentum from the collapse. However, 
with the other sign of f]H the magnetic braking is increased 
so that more angular momentum is transported to the en- 
velope. As magnetic braking due to Hall diffusion does not 
stop acting once no angular momentum remains (as am- 
bipolar diffusion does) it could also then spin the collapsing 
fluid back up in the opposite direction to the initial rota- 
tion. This acceleration is only possible with Hall diffusion, 
and it has the potential to completely resolve the magnetic 
braking catastrophe. 

Because of its tendency to move the magnetic fleld in 



© 2012 RAS, MNRAS 000.[Tlf24l 



18 C. Braiding and M. Wardle 



unusual directions Hall diffusion is usually overlooked in sim- 
ulations of gravitational collapse and star formation. It has 
been shown that the Hall effect is important to the dynamics 
of the collapse, particularly the magnetic braking behaviour 
which determines the existence and size of the rotationally- 
supported protostellar disc. The handedness of the response 
of the collapse to the inclusion of the Hall effect has obvious 
dynamical and potentially observable consequences for the 
gravitational collapse of molecular cloud cores, which must 
be studied more closely if the dynamics of the star formation 
process and the variations observed across YSOs and their 
discs are to be properly understood. 
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APPENDIX A: VERTICAL AVERAGING 

In order to produce a set of disc equations that depend only upon r and t, we average Equations [7 til2l vertically over the 
disc in order to reduce the dimensionality of the problem. As the equations for the conservation of mass, radial and angular 
momentu m and the vertical hydrostatic balance remain unchanged from KK02, the reader is directed to their appendix A 
or §2.3 of lBraidinid l|201lh for the details of this averaging, and present here only the derivations of the vertically-averaged 
induction equation and the azimuthal field component. 



Al z-component of the induction equation 

The z-component of the induction equation, given in Equation 1121 is expanded more completely as 



dt 



ld_ 

r dr 



B2 



dz dr 

2 , o2- 



yrB.+v(^~^] + '-§{B.^ + ^i-(rB 



B 



dB. 



Br d 



dz 



dB, _ dBr 
dr dz 



{Bt + B. 

The magnetic flux enclosed within a radius r is given by 

*(r-) = *c + 27r / B4r')r'dr' 



dr 



dB^ B^B^ d , , 

BrB^— 1 — {rB^^ 

oz r or 



where ^'c is the fiux within the central point mass. This equation is then rewritten in differential form as 

1 a* 



B, = 



and its derivative with respect to time is 



dB, 
dt 



27rr dr ' 

1 d /a* 

27rr dr \ dt 



(Al) 



(A2) 



(A3) 



(A4) 



This is substituted into Equation I All and the partial derivative with respect to r and the factor of r ^ are cancelled to obtain 

r VrB, + ri I 



1 a* 
2^~dt 



dBr _ dB, 
dz dr 



Vh 
B 



dB, Br d 



rjA 

£2 



{B^ + Br 



dB, _ dBr 
dr dz 



dB^ B^B^ d , . 

Br B4, — \ — {rB^ ) 

dz r dr 



(A5) 



The T] and riH,A terms depend on B'''^'^ respectively, so the leading fractions of the diffusive terms may be ignored as the 
integration over z is performed. 

The flux, magnetic force and Ohmic diffusion terms are integrated over the disc height to give: 



f 

J — C 



1 a^ 
2^ aF 



dz — 



_i_a* 

2^~dt' 



+H 



2H a^ 
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(A6) 



VrBzdz = 



VrB^Z 



+ H 



= 2HVrB, 



and 



+ °° ^ dBr 

dz 



dB, 
dr 



dz ■ 



-H dz 
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dr 



dz = 2{ B, 



H 



dB, 
dr 



The Hall diffusion terms are rearranged into the form 



^^^_^B^a a 

dz r dr dz 



r, dB Br d 

B^— 1 — (r_B0 ; 

dz r dr 



(A7) 



(A8) 



(A9) 



and the vertical scaling of the azimuthal field component is substituted into the first term and the solenoidal condition 
(Equation [TT} is applied to the second: 



dBi 



Br a 



^ dz ~^ r dr^ dz 



a / BzB^^sz\ B^ a 



H 
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The integral of the Hall terms may then by written as 

dz — 



f 

J — c 



■ dB^ Br d 

Bz—^ 1 TT-yrB^ 

dz r dr 
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al V H 
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dz, (All) 
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which, after the vertical scalings of Br and are substituted into it, is evaluated to give: 



dB^ 
dz 



B^d_ 
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B,B. 



z-^ d fr^Br.sB. 
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2H^ d fr'^BrsBA, 



3r2 dr 



(A12) 



Finally, the ambipolar diffusion terms are expanded out and integrated. The first of these terms is straightforward, as 
B^ is regarded as constant with height unless specifically differentiated with respect to z and may be taken outside of the 
integral, which is solved to obtain 



{Bl + Bt, 



dr 



dz • 



dB, 
dr 



+H 



Bi\dz = 2H 



dB, [Bl 
dr 



(A13) 



The second of the ambipolar diffusion terms is rearranged into the form 

dB dB d dB 

{Br + B^)— = Br—. h Bz-^{BrBz) — BrBz— — 

oz oz az dz 

to which the solenoidal condition (Equation Iff |l and the scalings for the other field components are applied. The integral of 
this term is then 



(A14) 
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(A15) 



this is evaluated over the height of the disc to give 

?3 'A 



/ + 00 
{Bl + Bl 
- oo 



dBr 

dz 



dz 



+ 



H 



+ 



3rH dr 

" d 



2 3 2 2 2 2 

gBr.s + 2B^Br,3 + -H Br^sB, 



dr 



H 

[\n{rBr 



dr 



[InH] 



(A16) 



The third of the ambipolar diffusion terms is again straightforward; it is vertically-averaged by applying the vertical scalings 
to the radial and azimuthal components to the field and then performing the integral over z to find 



' BrB^^dz^ 
dz 



dz 



H 



2 2 
dz — -Br,aB^^, 



(A17) 



Finally, the last of the ambipolar diffusion terms in Equation IA5I is averaged by substituting in the vertical scalings of the 
field components and then performing the integral: 
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Collecting all of these integrated terms into the same order as in Equation IA5I then gives the full vertically-averaged 
induction equation: 
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It is clear from this equation that the azimuthal field is pivotal in causing Hall drift in the radial direction; B,^,s should not 
be neglected, even in axisymmetric models. 



A2 Azimuthal field component 

The vertical angular momentum transport above and within the pseudodisc is achieved by magnetic braking, especially during 
the dynamic collapse phase inwards of the magnetic diffusion shock. It is assumed that magnetic braking remains the dominant 
angular momentum transport mechanism during the subsequent evolution of the core, although it is likely that a centrifugally- 
driven disc wind may domi nate in the innermost Keplerian disc. The approach to modelling the magnetic braking adopted 
here is adapted from that of lBasu &: Mouschoviaj ([1994) for the pre-point mass formation collapse phase. This formulation is 
not well-defined in the innermost rotationally-supported regions of the disc, where the calculated magnetic braking becomes 
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stronger than is expected and the angular momentum transport is expected to be dominated by a disc wind (this is discussed 
in more detail in !j5]). A cap is then placed upon the azimuthal magnetic field component in order to ensure that it does 
not greatly exceed the vertical component; because of this the magnetic braking prescription is not expected to introduce 
significant errors into the inner regions of the calculations. 

External to the pseudodisc the magnetic field is considered to be frozen into the low-density, constant-pressure external 
medium, which has density pcxt and angular velocity Q^. Within the external medium the magnetic field assumes the value 
B = ByciZ, and the exterior fiux tubes corotate with the core. Because the transition region has a low moment of inertia 
relative to the core, and the crossing time for Alfven waves is always much sm aller than the evolutionary time of the core, 
the transition region can relax to a steady state during all stages of contraction (|Basu fc Mouschovias|[l993 ). 

The induction equation under IMHD implies 

{Bp ■ V)Q = (A20) 

where Bp is the poloidal field, so that the angular velocity Q is constant on a magnetic surface. The force equation is similarly 

{Bp ■ V)rB^ = 0, (A21) 

which further implies that rBtf, does not change along the field lines. The neutral particles carry the torque and angular 
momentum is carried upwards by torsional Alfven waves generated by the rotation of the disc. 

Over a period of time dt an amount of material equal to 2-n prrcfdrrcf moves from the undisturbed position r-rddr-rcf in the 
external medium along a flux tube with angular velocity to a radius rdr at the disc surface. The angular momentum of the 

gas goes as 

dj = -[27r/9oxtr-rcfdrrcf](VA,oxtdt)rrcf(^^ - ^b), (A22) 

where VA.cxt, the external Alfven speed, is given by 

yA,ext = . (A23) 

For purely azimuthal motions in the external medium, the total angular momentum in each flux tube is conserved. This 
angular momentum must be removed from the disc at a rate equal to 

^ = —2TTr^cfVA,extPc^t{^ — flbWiafdrrci, (A24) 

dt 

_ 27rr^ef(yA,cxtPcxt)(^^ - 0,b)rre{drrei (A25) 

Tvrdr 

The amount of fiux remains constant along fiux tubes, so that the flux through the disc inside of a radius r is equal to the 
amount of fiux through the cylindrical external cloud inside of the radius ryct- 

/r 
2nr'B,,eq{r')dr' = TTr^,tB,,i, (A26) 



which gives a torque on the disc 



= 2-KrB^^e.qdr = 2-nrrcfBrcfdrrcf (A27) 

and 

fiefdVYei Bz^eq 



rdr B„ 



(A28) 



so that the torque in Equation I A25I becomes 

2r^ 



^ ^, - nb)Bz,,q / B,,ip,^t \ ^ _ (n - Qt)Bz,eq{^/2-K) ^^^Q) 



The torque per unit area on the disc is given by 

N = !:^f£2^; (A30) 
2tv 

combining Equations IA29I and IA30I gives the steady state azimuthal magnetic field component at the surface of the disc: 

(equation 26 of Basu & Mouschovias, 1994; equation 3 of KK02). It is clear that the properties of the external medium 
determine the conditions at the disc surface. This steady state approximation requires that the ratio of the Alfven travel 
time in the external medium to the initial radius of the cloud be less than the evolutionary timescale, which scales with r as 
■~ r/|l/r|. For the rotationally-supported discs presented here \Vr\ < Cs (and — > as r — > 0), which is much smaller than 
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the Alfven speed for the adopted gas temperature of lOK (VA.ext ~ 5cs); this imphes that the assumption of rapid braking of 
the core should not introduce large errors into the solutions. 
The angular velocity Q, is given by the equation 



where, using r]p = r]A + rj, 



n^-iV4> + VB4,), 

r 



Vb4. = \nH{V X B)x - ??p(V X B) 



X B 



(A32) 



(A33) 



This equation is then expanded out to become 



B2 



^ ( {Bi + B,^) ( ^ - ^ ) + B,B.-^i-(rB. 



B 



dB, dB 



VP 

B2 



dr 



{Bl + Bl + Bl) 



dz 

Br_d_ 

r dr 



1 d 



[rB^,) + B 



r dr 

dB^ 
dz 



Bd, B, 



dB^ 
dz 



(A34) 



Most of the terms in this equation have direct analogies in Equation IA5I and the individual steps of the vertical integration 
are not reproduced here. The vertical averaging gives 
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This equation is simplified as in 
final form of Vb4, is then: 



I by omitting any terms of order 0{H/r) save for the [Br,s — H{dBz/dr)] term; and the 
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This is equivalent to the ion-neutral drift velocity adopted by KK02 (their equation 9), with the inclusion of terms describing 
the effect of Hall diffusion. 

The ilfc term is dropped from Equation IA31I as the molecular cloud rotation rate is slow compared with that of the 
collapsing material. Rotation is dynamically important in the inner regions of the solutions presented in this thesis, while it 
is not important in most molecular clouds, so it is reasonable to declare that Q ^ fit and dismiss Qb as small. The external 
Alfven speed, VA,ext, is treated as a constant with respect to the isothermal sound speed in these calc ulations, param eterised by 
the constant a (defined in Equation [T^} . This scaling of VA,cxt is reasonable as the observations bv ICrutcheJ (| 199il ) i 



that Va ~ 1 km s ^ over at least four orders of magnitude in density (~ 10^-10^ cm 
Equations IA36I and If 91 are substituted into [TS] to find that 



ndicated 

in their observed molecular clouds. 
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Note that B has an implied B^^s dependence; this is typically solved for numerically when calculating the azimuthal field. 

For the inner solutions, increases with decreasing r (proportional to r"^''^); this would make B^^s the dominant field 
component at the surface near to the central point mass. Such behaviour is not expected in a real disc, where internal kinks 
of the field and magnetohydrodynamical instabilities (for example, the magnetorotational instability) should reduce the value 
of B^ at the surface. An artificial limit on B^.s is imposed: 
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(A38) 



where 5 is a parameter of the model usually chosen to be (5 = 1 in order to ensure that the azimuthal field component 
does not exceed the vertical component. KK02 point out that this value quite conveniently corresponds to that expected for 
a rotationally-supported disc where the vertical angular momentum transport is dominated by a centrifugally-driven wind. 
Applying this cap to Equation I A37I then gives the final equation for B^^j,: 
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(A39) 
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